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Molecular dynamics simulations of Lennard- Jones systems are performed to study the effects of 
dissolved gas on liquid-wall and liquid-gas interfaces. Gas enrichment at walls is observed which for 
hydrophobic walls can exceed more than two orders of magnitude when compared to the gas density 
in the bulk liquid. As a consequence, the liquid structure close to the wall is considerably modified, 
leading to an enhanced wall slip. At liquid-gas interfaces gas enrichment is found which reduces the 
surface tension. 

PACS numbers: 68.08.-p,68.03.-g,68.15.+e,83.50.Rp 
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The precise determination of the hydrodynamic 
boundary condition, slip vs. no-slip, is currently a mat- 
ter of active debate. A growing number of studies, ex- 
periments 0, Q, as well as simulations 0, 0, 0, 0, @ , 
strongly indicate that the classical no-slip condition, a 
more than 200 year old dogma, is violated. Though it is 
difficult to identify clear trends, the investigations sug- 
gest that increasing hydrophobicity and an increasing 
amount of dissolved gas in the liquid favor larger slip. 
Note, however, that slip has been reported for hydrophilic 
surfaces as well Q. 

Despite many investigations, slippage behavior and 
its origin are far from being understood. A possible 
cause is the presence of so-called surface nanobub- 

bles, i.e., nanoscale bubbles located on a solid surface 
that is immersed in liquid. Many recent experiments 
support the notion of surface nanobubbles, in particu- 
lar atomic force microscopy measurements [lfjj, but also 
other techniques 0]. Similar to the trends for wall slip, 
hydrophobicity and dissolved gas favor nanobubbles. For 
gas-saturated liquid nanobubbles are found on hydropho- 
bic surfaces, whereas usually nanobubbles are not ob- 
served for hydrophilic and/or degassed liquid, suggesting 
gas- rather than vapor bubbles. In spite of growing ex- 
perimental evidence for their existence, it is unclear how 
and why they form and why they are apparently stable. 

Other examples for the interplay between hydrophobic 
interfaces and dissolved gases are colloidal suspensions 
and emulsions where the stability is considerably 
influenced by the presence of dissolved gases. Moreover, 
recent neutron reflectivity measurements reveal a de- 
pendence of the width of the hydrophobic wall-water in- 
terface on the amount and type of dissolved gas. 

Though the above mentioned experiments clearly 
demonstrate the importance of dissolved gases for the hy- 
drophobic wall-liquid interface, a profound understand- 
ing is still lacking. Molecular dynamics simulations are 
a promising approach to address this issue. However, 
previous simulations, for instance of slippage 0,0,00]; 
were restricted to pure liquids without dissolved gases. 
How do gases effect liquid-wall interfaces? How do the 
effects change with hydrophobicity or for different gases? 



Is wall slip enhanced? It is the aim of this Letter to ad- 
dress these issues by means of molecular dynamics simu- 
lations. Control parameters are the amount of dissolved 
gas, the hydrophobicity of the wall, and the type of gas. 
Liquid-gas interfaces, which serve as reference systems 
and are important in their own right are studied as 
well. 

Simulations are performed for fixed particle number, 
volume and temperature T=300 K using the GROMACS 
code Periodic boundary conditions (p.b.c.) are ap- 
plied in x, y— and z-direction. Three different particle 
species (liquid/gas/wall) with mass m=20amu are sim- 
ulated. Liquid and wall particles have the same molecu- 
lar diameter <7=0.34nm. Particles interact via Lennard- 
Jones (LJ) 6-12 potentials with a cutoff r c =5a, which is 
larger than the value 2.5a usually applied for bulk liquids, 
in order to account for inhomogeneities at interfaces. The 
energy scale eu for liquid-liquid interactions is fixed to 
eiiK,\.2k,BT with Boltzmann's constant fcs. To model 
an inert gas without wall affinity the energy scales for 
gas-gas and gas-wall interactions are e gg —e gw ^0AkBT, 
which is close to e gg of Argon. The temperature T is be- 
low (above) the critical temperature T c of the liquid (gas) 
particles 0,0. The time step is di=0.005r with the 
characteristic time T=cr-\/ra/e;;«0.9 ps. During produc- 
tion runs, the simulations are weakly coupled to a heat 
bath using the Berendsen thermostat ^g] with a relax- 
ation time tt=10t. A perfectly stiff wall is simulated by 
solid particles that are frozen on a fcc-lattice with den- 
sity p w ~0.96a~ 3 . The center of mass velocity is removed, 
apart from the flow simulations. 

Four microscopic control parameters (i)-(iv) are tuned 
which change the properties of interest. To simulate dif- 
ferent gases (i) the energy scale e g i for gas-liquid inter- 
actions is varied, as well as (ii) the molecular diameter 
a g of the gas particles (a g i=a gw =0.5(a g +a) is applied). 
Expressing e g i and a g in terms of fc^T and a, respec- 
tively, the combinations (e g i,a g )—(0A,l), (0.692,1.47), 
(0.692, 1.62), and (0.712, 1.62) are studied, which are de- 
noted as gas types (A) — (D). To identify effects due to 
the gas (iii) the number of gas particles N g is changed 
from N g =0 (pure liquid) to the finite value N g —22%. The 
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FIG. 1: (color online) Starting configurations to study (a) 
liquid-gas interfaces, and (b,c) liquid films at walls. Initially 
the particles are located on a lattice, which 'melts' during 
equilibration, forming a liquid while the vacuum is filled by a 
vapor-gas phase. Due to p.b.c. in (b,c) the wall terminates the 
vapor-gas phase in z-direction. The scale is given by Z«16<r. 



hydrophobicity of the wall is varied by (iv) the ratio 
fito/cH with the scale ei w for liquid-wall interactions, en- 
abling simulations of hydrophilic and hydrophobic walls. 
These microscopic parameters (i)-(iv) determine macro- 
scopic properties such as gas solubility, gas concentration, 
surface tension, and the contact angle. 

Before addressing liquid-gas mixtures at walls it is 
worth discussing liquid-gas mixtures without walls. Ini- 
tially, liquid and gas particles are located on a lat- 
tice ('fluid cube') in the center of a rectangular sim- 
ulation box, Fig. H[a). After an equilibration period 
(9xl0 6 (ii«40ns), which consists of a series of subsequent 
microcanonical simulations at T=300 K, the system is in 
a steady state with a liquid film perpendicular to the 
z-axis, in coexistence with the vapor-gas phase. The to- 
tal number of particles is iV=2916. Fig. presents den- 
sity profiles obtained from time averaging (10 6 dt after 
equilibration). One can clearly observe an enrichment 
of gas in the interfacial region, before the gas density 
falls off towards its value in the bulk liquid (similar ob- 
servations have been made for liquid-liquid mixtures [lflj ] 
with a much stronger attraction between different parti- 
cle species). A gas particle close to the interface expe- 
riences attractive forces from particles in the vapor-gas 
phase as well as in the liquid film. Since the density 
in the liquid film is much larger than in the vapor-gas 
phase, the resulting force is directed towards the liquid 
film, which leads to the nonmonotonous density profiles, 
even for gases with small gas-liquid interactions as for 
(A). Note that the amount of gas in the bulk of the liq- 
uid is similar for all gases, (6.25±2.75) x10~ 4 ct~ 3 , though 
the importance of different factors involved in the pro- 
cess of gas solution are expected to differ for the gases 
(A) — (D). To illustrate this, consider the energy scale e 9 ; 
of (D) which will facilitate gas solution compared to (A) 
with its smaller value of e s ;, but this is counteracted by 
the larger size of the (D) particles. 

Does gas change the surface tension 7? Experiments 
show that gases decrease 7 (which has been proposed 
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FIG. 2: (color online) Density profiles for liquid-gas mixtures 
at phase coexistence. Gas enrichment at the interface can 
be observed (see inset). The liquid profiles were similar for 
the different gases, wherefore only one is displayed. Varying 
parameters for the gases (A) — (D) are (e g i/e gg ,a g /a) = (l, 1), 
(1.73, 1.47), (1.73, 1.62), and (1.78, 1.62). 



to be crucial for bubble nucleation), but the reason 
was stated to be unknown |14J. Applying the standard 
Kirkwood-Buff 20] formula to calculate 7, a decrease 
of 7 due to the gas adsorption is found as well. Ex- 
pressed in terms of e;;/cr 2 , the average value of 7 for the 
liquid-vapor interface (N g =0) is 7«0.74 (in agreement 
to typical values for LJ-fluids 0]), which is reduced to 
approximately (0.72,0.6,0.56,0.56) for (A) -(D). Addi- 
tional simulations show that the decrease of 7 is enhanced 
for increasing gas pressure, just as in experiments. 

What changes in the presence of walls? Since the con- 
tact angle 9 is of vital importance, the walls are first 
characterized by simulations of droplets at walls. There- 
fore, a 'fluid cube' composed of liquid particles is initi- 
ated on a wall. After dynamical evolution (2.5xl0 6 dt) 
the density profiles (obtained by 10 6 dt time averaging) 
allow to estimate 9 numerically, see Fig. [3J The trend of 
the obtained contact angle with the hydrophobicity pa- 
rameter ei w /eii is consistent with what one would obtain 
from the rough estimate cos#l~— 1+2 (p w ei w )/(pi e«) 
(with densities p w and pi of wall and liquid) based on the 
Laplace expression of surface energies |20j . 

What is the molecular structure of liquids in contact 
with walls, in particular in the presence of dissolved gas? 
With well controled walls in place, we proceed to investi- 
gate this issue. Therefore a 'fluid cube' of liquid and gas 
particles is initiated close to a wall, Fig. IHb) . The effect 
of hydrophobicity is studied by changing ei w /eu, as dis- 
cussed above. In order to probe the effect of dissolved gas 
we compare simulations with N g =0 (pure liquid) to sim- 
ulations with N g =228 for the gases (A) — (D). The num- 
ber of liquid particles is Ni— 2688. After an equilibrating 
period (9xl0 6 di and 12AxlO e dt for the hydrophilic and 
hydrophobic wall) the system consists of a fluid film in 
phase coexistence on one side and which is in contact 
with a wall on the other side. The left part of Fig. 0] 
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FIG. 3: (color online) Liquid density profiles of droplets at 
walls, to characterize the walls. Tuning the microscopic at- 
traction ratio ei w /eu (with fixed en) results in a change of 
the macroscopic observable 6. Left, (a): hydrophilic wall, 
eiti>/e«=0.617 with a measured contact angle 9~80° . Right, 
(b): hydrophobic wall, ei m /eu=0.333 with 0«115°. The x- 
extension of the initial droplets ('fluid cube') equals the x- 
extension of the simulation box, leading to (hemicylindrical) 
droplets, translationally invariant in x-direction. 



shows liquid and gas density profiles (10 6 dt time averag- 
ing) close to the hydrophilic wall. There is a tiny increase 
of the gas density at the wall, which is too small to be 
observed on the scale of Fig. 0] The pure liquid exhibits 
the usual layering [5J, Il7l l21| close to the wall, which is 
hardly altered by the presence of the gas. 

How do the density profiles change for a hydrophobic 
instead of a hydrophilic wall? Fig. 0^b) shows liquid and 
gas density profiles close to the hydrophobic wall. One 
immediately observes a dramatic increase of the gas den- 
sity in the vicinity of the wall. For gas (A) the density in- 
creases by a factor «50 when compared to the gas density 
in the bulk liquid, and the gas enrichment is even more 
than two orders of magnitude for the gases (B)~(D). 
Furthermore, the liquid structure close to the wall is 
drastically changed. The pure liquid exhibits layering, 
which is less pronounced than for the hydrophilic case, 
Fig. 01a). The liquid structure is only slightly altered by 
gas (A) but it is greatly diminished by the presence of 
gases (B) — (D). The gas enrichment leads to a consid- 
erable reduction of the liquid density in the vicinity of 
the hydrophobic wall. We stress that for the hydropho- 
bic wall for all gases e gw =ei w , which shows that the gas 
enrichment is not caused by a strong gas-wall interac- 
tion. The gas enrichment is associated with a reduced 
diffusion of the gas perpendicular to the wall. Since the 
simulations leading to Fig.0]started with a high gas con- 
centration in the liquid, Fig.^b), a proper equilibration 
is a delicate issue (for gases (B) — (D)). Therefore we con- 
firmed the results by simulations starting from a contrary 
initial configuration, Fig.^c). Here, the gas particles are 
initially completely separated from the wall by the liq- 
uid film. The gas enrichment obtained after 12.4xl0 6 dt 
equilibration from simulations initiated from Fig. ^c) is 
71%, 56%, and 92% of the gas enrichment for the gases 



(B)—(D) depicted in Fig. gb) (23. Hence, the tremen- 
dous gas enrichment as well as the considerable reduction 
of the liquid density at the wall are reproduced even in 
simulations starting from the configuration Fig.^c). 

The width of the region of gas enrichment is of the 
order of a, similar to the observations in 13]. Thus, we 
basically find a monolayer of gas particles adsorbed at 
the wall. Though the gas enrichment shows some re- 
semblance of surface nanobubbles it is still different from 
gas bubbles with heights of several nanometers as exper- 
imentally observed [10|. Interesting objectives for future 
research are to clarify if the gas enrichment constitutes a 
reservoir for nanobubbles and if the reduced diffusion of 
gas perpendicular to the wall helps stabilizing them. 

What causes the dramatic gas enrichment? Energeti- 
cally the system benefits from gas-liquid interactions due 
to € g i, but gas particles in the bulk liquid occupy space 
due to cr g which is unfavorable (reduction of liquid-liquid 
interactions). Gases at the liquid interface, however, re- 
duce the energy with little disturbance of the liquid-liquid 
interactions, and energetic contributions from ei w , which 
are diminished, are small for hydrophobic walls. Accord- 
ing to this explanation, the gas enrichment increases with 
increasing e g i and a g (for similar gas concentration in the 
bulk liquid), just as observed in the simulations. 

To probe the effect of the gas on the slippage behav- 
ior, we apply a constant force f y =2.27xlO~ 3 eu/<T in y- 
direction (parallel to the wall). As usual only the veloc- 
ity component perpendicular to the flow (x-component) 
is thermostated (Langevin thermostat with tt=t). 
The unforced systems after their equilibration phase (ini- 
tial configuration Fig.^c) for (B) — (D)) are further equi- 
librated (1.2xl0 6 cft) while applying f y , and hereafter 
production runs (10 e dt) yield velocity profiles shown in 
Fig. El The velocity in the liquid film at the hydrophilic 
wall is not altered by the gas. Contrarily, at the hy- 
drophobic wall the gas significantly changes the veloc- 
ity profiles, leading to an increase of the average ve- 
locity. Estimates of the slip length X=\v v /d z v v \ wa ii 
using the fits depicted in Fig. 03 (dashed lines) yield 
A«(3.7, 3.4, 4.5, 7.0, 7.9)a for the pure liquid and liquid in 
the presence of (A) — (D), respectively. Hence, the pres- 
ence of gas can significantly increase the slip length. 

In conclusion, our results support the experimental 
findings that gases dissolved in liquids, although present 
only in low concentration in the bulk liquid, can have a 
strong influence on the structure of the liquid-wall inter- 
face, due to gas enrichment at hydrophobic walls. Future 
studies of phenomena associated with the hydrophobic 
wall-liquid interface therefore must take dissolved gases 
into account. This holds, e.g., for the appearance of 
nanoscale bubbles^the study of slippage, and the break- 
age of nanofilms [23|. 
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FIG. 4: (color online) Liquid and gas density profiles (same symbols for liquid and gas) for liquid films in contact with walls. 
Left, (a): hydrophilic wall, the liquid exhibits the usual layering which is not altered by the the presence of the gases. Right, 
(b): hydrophobic wall, note the tremendous increase of the gas density and the greatly diminished liquid density at the wall. 
Varying parameters for the gases (A) -(D) are (e gl /e gg , cr g /cr)=(l, 1), (1.73, 1.47), (1.73, 1.62), and (1.78, 1.62). 
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FIG. 5: (color) Velocity profiles of the liquid close to the wall. 
The gas enrichment at the hydrophobic wall clearly enhances 
the average velocity, contrary to the hydrophilic case, where 
the gas has no influence. The dashed lines are quadratic fits 
which are used to estimate the slip length (see text). 



